Multiband-Driven Superfluid-Insulator Transition of Fermionic Atoms in Optical 
Lattices: A Dynamical Mean-Field-Theory Study 
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The superfluid-insulator transitions of the fermionic atoms in optical lattices are investigated by 
the two-site dynamical mean-field theory. It is shown that the Mott transition occurs as a result of 
the multiband effects. The quasiparticle weight in the superfluid state decreases significantly, as the 
system approaches the Mott transition point. By changing the interaction and the orbital splitting, 
we obtain the phase diagram at half filling. The numerical results are discussed in comparison with 
the effective boson model. 
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Since the superfluid of trapped atomic Fermi gases 40 K 
and 6 Li was observed 0, 2, El, 0, 0, @| , intense theoretical 
and experimental studies have been done for ultracold 
atomic Fermi gases. In these experiments, magnetic-field 
Feshbach resonances provide the means for controlling 
both the strength of the interaction between fermionic 
atoms and its sign Q • These tunable interactions enable 
us to observe the crossover between the BCS superfluid 
for the weak attractive interaction and the Bose-Einstein 
condensate (BEC) of bound pairs for the strong attrac- 
tive interaction. Furthermore, by loading atoms into op- 
tical lattices, diverse interaction configurations can be 
introduced. The combination of these two experimental 
techniques plays an important role in the study of ultra- 
cold fermionic atoms and offers the experimental descrip- 
tion of various intriguing quantum many-body phenom- 
ena. 

Recent experiments revealed fascinating aspects of 
fermionic atoms in three dimensional optical lattices 
[H, In the ETH experiment, a band insulator in the 
lowest band was produced, i.e. two atoms in different 
hyperfine states occupy the lowest state per lattice site 
[8|. Controlling the interaction, they further observed 
the partially populated higher bands. In the MIT ex- 
periment, a superfluidity of fermionic atom pairs in an 
optical lattice was observed [i| . By increasing the depth 
of the lattice potential near the Feshbach resonance, a 
superfluid-insulator transition was observed. They ar- 
gued that the insulating state was the Mott insulator. In 
these experiments, it was argued that the usual single- 
band Hubbard model was no longer applicable, because 
the strength of the on-site interaction exceeded the gap 
between the lowest and the next-lowest bands. For de- 
tailed investigations of the experiments, accordingly, the 
effects of the higher bands have to be taken into account. 

Stimulated by these experiments, several theoretical 
studies on the superfluid-insulator transition were carried 
out 13, II , 12j • However, both the correlation effects and 
the multiband effects have not yet been investigated well 
enough to discuss the transition from superfluid to Mott 
insulator. Precise studies including both effects are thus 



required. 

In this paper, we investigate the superfluid-insulator 
transition of interacting fermionic atoms in optical lat- 
tices, taking into account the multiband effects. For this 
purpose, we make use of a dynamical mean-field theory 
(DMFT) 0. This method enables us to treat local cor- 
relation effects precisely. We show that the multiband 
effects cause the transition from superfluid to Mott in- 
sulator. The phase diagram over the interaction from 
repulsive to attractive at half filling is presented. Fea- 
tures of each phase are discussed. 

We consider that the fermionic atoms are in a periodic 
optical lattice potential. In the low-tunneling limit, each 
lattice potential is regarded as a harmonic one [HQ. In 
this case, the lowest orbital is nondegenerate, while the 
next-lowest orbital is three-fold degenerate. We approx- 
imately model the system as follows. The degeneracy 
of the next-lowest orbital is neglected for simplicity. We 
introduce the effects of the lowest and the next-lowest or- 
bitals. The hopping integrals between the same orbitals 
of the nearest neighbor sites are assumed to take the same 
value for both orbitals. According to the usual experi- 
ments, we assume that the system involves the same num- 
ber of fermionic atoms in two different hyperfine states. 
These two hyperfine states are described as the pseu- 
dospin. The simplified model Hamiltonian we consider 
here reads 

H = (t - vSi,j)4 a „c jaa + —^2(n i2t7 - n ila ) 
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where Ci aa is the fermionic annihilation operator for the 
state with pseudospin cr(=f ,|) on orbital a(—l, 2) in the 
ith lattice site and n iacr — c[ aiy c iaa . t represents the 
hopping integral, fj, the chemical potential, D the split- 
ting between the two orbitals. For the lowest and the 
next-lowest orbitals, the intraorbital on-site interaction 
U between fermionic atoms, and the interorbital ones U' 
and J fulfill the relation U' = J = U/2. The Hamiltonian 
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is known as a two-band Hubbard model. We assume that 
the intraorbital attractive interaction induces an s-wave 
superfluid state. 

We investigate the ground state by using DMFT. In 
DMFT, the lattice model is mapped onto a single im- 
purity model connected dynamically to a heat bath. We 
solve it self-consistently [13j . This retains nontrivial local 
quantum fluctuations missing in conventional mean-field 
theories. There are various methods to solve the effective 
impurity model. In fact, the single-band attractive Hub- 
bard model has been analyzed by some DMFT methods 
3 EE El El, HI. We apply here the two-site DMFT 
method [21(, which allows us to discuss the Mott transi- 
tions of orbitally degenerate lattice fermions qualitatively 



22l . 23f . Comprehensive investigation can be done, e.g., 
for the determination of phase diagrams. To study the 
superfluid of lattice fermions, we extend this method to 
the case when the superfluid order exists. The exten- 
sion of the two-site DMFT to the multiband system is 
straightforward. 

We use a semicircular density of states (DOS) p(uj) = 
4/(itW)^1-A(lu/W) 2 , where W is the band width. 
Since the hopping integrals are assumed to be indepen- 
dent of a, W — At and the DOS are the same for both 
bands. The chemical potential is set /i = U/2 + U' — J/2 
so that particle-hole symmetry can be satisfied in both 
bands. In this case, the system is half filling. We neglect 
the density wave state which degenerates with the super- 
fluid state in the symmetric case, because we focus on the 
superfluid state. In the following, the hopping integral t 
is used in units of energy. 

We calculate the quasiparticle weight Z, the superfluid 
order parameter $ = {c ia ic ia ^), and the filling of each 
band ni, rii for U' = J = U/2 in the attractive U. Z 
represents the coherent spectral weight of the Bogoliubov 
quasiparticle [l8| , which enables us to examine the Mott 
transition accurately. For half filling, <!> becomes inde- 
pendent of i and a in the DMFT procedure. In Fig. 
QJa), we show the results for D = 1.0. As \U\ increases, 
Z decreases to a small value close to zero at U = —3.0. 
$ first increases and then decreases in \U\ > 2.7. At 
U = —3.0, $ vanishes discontinuously. The significant 
decrease of Z in the superfluid state is a manifestation of 
the strongly renormalized Bogoliubov quasiparticle. Be- 
cause of the finite D, n\ is larger than in \U\ < 3.0. 
As \U\ increases, m decreases and n 2 increases, and both 
change to 1 suddenly at U = —3.0. The results indicate 
that the quantum phase transition from the superfluid to 
the Mott insulator takes place at U = —3.0 for D = 1.0. 
The changes of the above quantities around U = —3.0 are 
rather discontinuous, suggesting that the transition may 
be of the first order. We also show the DOS of the low- 
est band for several values of U. We find a spectral gap 
around u) = 0, which is caused by the superfluid order. 
As \U\ increases, the spectral gap grows and the weight of 
the coherent part is reduced. The results are consistent 
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FIG. 1: The quasiparticle weight Z, the superfluid order pa- 
rameter and the filling of each band 711,712 for U' — J = 
U/2. The DOS of the lowest band are shown in the right pan- 
els of (a) and (b). The results show three types of transitions: 
(a) the transition between the superfluid and the Mott insu- 
lator, (b) the transition between the superfluid and the band 
insulator, and (c) the transition between the band insulator 
and the Mott insulator. 



with the decrease in the quasiparticle weight Z. Further- 
more, the weight of the incoherent part increases steadily. 
At U = —3.0, the coherent part vanishes and the spectral 
gap by the superfluid order around u = turns into the 
Hubbard gap. The results are consistent with each other 
concerning the quantum phase transition. Accordingly, 
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we conclude that the transition from the superfluid to 
the Mott insulator takes place at U = —3.0. We empha- 
size that the Mott transition is caused by the multiband 
effects. In fact, the transition from the superfluid to the 
Mott insulator has not been obtained in the single-band 
attractive Hubbard model ll| H, 2pj . 

In Fig. [IJb), we show Z, $, n%, and n-i as a func- 
tion of D for U = —1.5. The results are in noticeable 
contrast with those shown in Fig. [TJa). As D increases, 
Z increases to 1 and $ decreases to at D = 3.2. We 
also observe that n\ and ri2 change linearly from I in 
small D, and go to 2 and at D = 3.2, respectively. 
These changes are rather continuous, suggesting that the 
transition is of the second order. The results show that 
the quantum phase transition from the superfluid to the 
band insulator occurs at D = 3.2 for U = —1.5. This 
conclusion is confirmed by the DOS of the lowest band. 
As D increases, the spectral gap by the superfluid order 
becomes narrower and the weight of the incoherent part 
get smaller. For D — 3.2, the spectral gap vanishes and 
the band insulator emerges. 

In Fig. [He), we show the results for Z and m, in 
relatively strong attractive interaction U = —4.5. As D 
increases, Z changes from small value to 1, and the fillings 
change from ni = n-i = 1 to m = 2 and n% = 0. The 
results imply that orbital fluctuations are suppressed and 
the direct transition from the Mott insulator to the band 
insulator takes place. In these quantities we observe the 
hysteresis, which is a proof of the first order transition. 
Judging from the present results and those in Fig. Ufa), 
we conclude that the Mott transition for the attractive 
interaction is of the first order. The transition point is 
determined at D — 2.3 for U = —4.5, where the lowest 
energy is provided. 
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FIG. 2: Phase diagram for U' = J = U/2 at half filling. In 
the attractive side the bosonic fermion pair occupies either of 
the two orbitals in each site, while in the repulsive side the 
fermionic atoms occupy the respective two orbitals in each 
site. 



Changing D and U systematically, we investigate the 



quantum phase transition over the attractive and repul- 
sive U. For the repulsive U, we obtain [Z, 711,7x2) = 
(1, 2, 0) and (0, 1, 1). The former represents the band in- 
sulator and the latter the Mott insulator. We further 
obtain < Z < I and < n-z < I < n\ < 2, which repre- 
sent the metallic state. The results for U' = J = U/2 at 
half filling arc summarized in the phase diagram shown in 
Fig. [21 The band insulating phases smoothly connect be- 
tween the repulsive and attractive sides. The transition 
line of the band insulator agrees well with the result by 
the conventional mean-field theory: D = U/2 + 4. Com- 
bining the results discussed in Fig. [ljb), we conclude 
that the band-insulator transition is of the second order. 
As the attractive interaction becomes strong, the super- 
fluid state in the small D region (D < 2) changes directly 
into the Mott insulator owing to the strong renormaliza- 
tion of the correlation effects. In the appropriately large 
D region (2 < D < 4), by contrast, the superfluid state 
first changes into the band insulator which is suitable for 
the fermion pair, and then the band insulator changes 
into the Mott insulator which is favorable in large \U\. 

We now investigate the origin of the transition from 
the superfluid to the Mott insulator in another point of 
view. To this end, we map the model Hamiltonian at half 
filling to the effective boson model which is appropriate 
for the strong attractive U [261 ] . 



T^Boson — tcS 

<i,j> 



&& + 0&I>(ni-l), (2) 



where bi annihilates a bosonic fermion pair on the ith 
lattice site, = b\bi, tf s (= —2t 2 /U) represents the ef- 
fective hopping integral of the fermion pair, and U^ s (— 
—2U' + J — D) denotes the effective interaction between 
two bosonic fermion pairs in the same lattice site. In 
deriving the effective boson model, we assume that D 
is much smaller than \U\. Since the original model has 
only two orbitals, at most two bosonic fermion pairs oc- 
cupy the same lattice site in the effective boson model. 
Under the condition for the lowest two orbitals: U' = 
J = U/2(< 0), we obtain U^ > 0. The repulsive 
U® s drives the superfluid state to the Mott insulator of 
bosonic fermion pairs. On the basis of the effective bo- 
son model, we discuss the Mott insulating states in the 
repulsive and attractive sides. Although the numerical 
results show no difference between them, the situations 
are different. For the Mott insulating state of the repul- 
sive side, it is considered from the numerical results that 
the fermionic atoms in each site occupy respective two 
orbitals J which is consistent with the results obtained so 
far 24j, [25j . For the Mott insulating state of the attrac- 
tive side, by contrast, the bosonic fermion pairs occupy 
either of the two orbitals in each site. Therefore, the av- 
eraged fillings take the same value rt\ — = I, which 
is seemingly the same results as the repulsive side. As 
D increases, orbital fluctuations in both sides play differ- 
ent roles. In the repulsive side, the metallic state sand- 
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wiched between the two insulating phases is stabilized 
by orbital fluctuations 24], |25| . In the strong attractive 
region, on the other hand, orbital fluctuations of fermion 
pairs are suppressed as mentioned before, leading to the 
direct transition between the two insulating states. 



We thank A. Yamamoto and M. Yamashita for valu- 
able discussions. KI was supported by the Japan Society 
for the Promotion of Science. Numerical computations 
were carried out at the Supercomputer Center, the Insti- 
tute for Solid State Physics, University of Tokyo. 




-1 .2 -3 -4 -5 -6 -7 -8 

U 



FIG. 3: The superfluid order parameter $ and the quasipar- 
ticle weight Z as a function of the attractive interaction U for 
several values of J in U' = U/2 and D = 0. 



The expression suggests further that for J = 2U' 
and D = the effective interaction vanishes as U^g = 
and the Mott transition disappears. To confirm this sug- 
gestion, we calculate the quasiparticle weight Z and the 
superfluid order parameter 4>, releasing J from the re- 
striction U' = J = U/2. The results for several values 
of J in U' = U/2 and D = are shown in Fig. GO We 
find that Z ^ and $ ^ for J = 2U'. The results 
demonstrate that the Mott transition indeed disappears 
in agreement with the suggestion based on the effective 
boson model. On the other hand, the Mott transition oc- 
curs for J 7^ 2U'. The agreement between the numerical 
results and the effective boson model indicates that the 
results based on the two-site DMFT and the description 
for the Mott insulator are reliable. 

In summary, we have investigated the quantum phase 
transition of the fermionic atoms in optical lattices on 
the basis of a simplified model. Using two-site DMFT, 
we have shown that the transition from the superfluid to 
the Mott insulator at half filling is caused by the multi- 
band effects. As a precursor of the Mott transition the 
renormalizcd superfluid state appears, which is charac- 
teristic of the interacting multiband fermionic atoms in 
optical lattices. The phase diagram among the super- 
fluid, the Mott insulator, the band insulator, and the 
metallic state at half filling has been determined. 
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